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Abstract: In supergravity where the gravitino is the hghtest supersymmetric particle 
(LSP), the next-to-hghtest supersymmetric particle (NLSP) decays to the gravitino with 
a naturally long lifetime ( 10^ — 10*^ s). However, cosmological constraints favour charged 
sleptons with lifetimes below a year as the natural NLSP candidate. For this scenario we 
report a method to accurately determine the slepton lifetime and SUSY cross-section from 
observation of the decays of sleptons trapped in the material comprising the main detector 
(ATLAS, CMS). A measurement of the lifetime to 5% is possible after 3 years at nominal 
luminosity and running conditions. This method is sensitive to the cosmologically preferred 
stau lifetime of ~ 37 days and does not require the use of ancillary trapping volumes. 
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1. Introduction 

Most studies of supergravity models assume the lightest supersymmetric particle (LSP) 
is a standard model superpartner, such as a slepton or neutralino. However, there are 
theoretically plausible scenarios where the gravitino is the LSP lll§[||l[§]||0®il 



lOd Here the gravitino is a "superWIMP" with a weak scale mass of ~ 100 GeV and 



couplings suppressed by the reduced Planck scale (RPS), Mrps = -v/SttG/v, where Gat is 
the gravitational constant. If the gravitino is the LSP, the next-to-lightest supersymmetric 
particle (NLSP) decays to its standard model partner and a gravitino. The NLSP is a 
weak-scale particle decaying gravitationally and so has a natural lifetime of the order of: 
^RPs/^weak ~ lO^-lO'^s {tj\[lsp)- This extraordinarily long lifetime renders the decay of 
the NLSP particle undetectable in a typical collider detector. 

The superWIMP is a good candidate for dark matter and as such the gravitino LSP 
scenario is constrained by measurements of non-baryonic cold dark matter density. The 



NLSP decays will deposit electromagnetic [12| and hadronic [O] energy into the universe, 



which may upset the successful predictions of standard big bang nucleosynthesis (BBN). 
Also, it is possible these NLSP decays may distort the cosmic microwave background (CMB) 
from its observed Planckian spectrum. In addition, photons produced in NLSP decays are 
subject to bounds on the diffuse photon flux. However, it turns out the superWIMP 
scenario can satisfy all these constraints. 
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Neutralino NLSPs are highly disfavored [0] [|ll|. However, slepton and sneutrino NLSP 
scenarios have been analyzed and found to be safe |6|[0 with the result that the most 
natural NLSP candidates are sleptons, particularly the right-handed stau. Cosmological 
constraints exclude the upper range of t^lsp since, if the decays occur in a colder uni- 
verse, then the decay products would not be effectively thermalized. The CMB and BBN 
constraints, for typical thermal relic NLSP abundances, provide an upper bound on NLSP 
lifetimes that excludes tmlsp above roughly a year. 

Late NLSP decays may even resolve the leading BBN ^Li anomaly by reducing the 
predicted ''Li abundance to the low values favoured by observations ^ Q. The NLSP 
lifetime that best resolves the ^Li anomaly is 3x 10^ s [^, of the order of a month. An- 
other motivation follows from considerations of leptogenesis [ |14[ . In the gravitino LSP 
scenario, the gravitino does not decay. Thus, the reheat temperature is bounded only by 
the overclosure constraint on the gravitino density. For a gravitino mass Mq ~ 100 GeV, 
reheat temperatures as high as ~ lO^'^ GeV are allowed consistent with thermal leptoge- 
nesis Further connections between gravitino LSPs and leptogenesis are discussed 



in Ref. 117 



The collider implications of a gravitino LSP with a charged slepton NLSP with life- 
time under (but not much under) a year have been investigated |l^. In particular, the 
possibility of trapping sleptons in material (water tanks) placed just outside Large Hadron 
Collider (LHC) or International Linear Collider (ILC) detectors was investigated. In this 
approach the material in which the sleptons are trapped is moved to an underground loca- 
tion (underground resevoirs) in order that slepton decays may be observed in a relatively 
background-reduced environment. At the Large Hadron Collider (LHC), the investiga- 
tors found that, depending on the mass scale of supersymmetry and the luminosity of the 
machine, tens to thousands of sleptons may be trapped each year. This implies that per- 
cent level studies of sleptons may be performed with resulting insights into supergravity, 
supersymmetry breaking, dark matter and dark energy ||19|. 



In this paper, we report a method of measuring the slepton liftime in a typical LHC 
detector modelled on ATLAS, using the decays of sleptons trapped in the detector itself. No 
use is made of additional water traps or underground resevoirs. The cosmic ray background 
is essentially eliminated by the depth of the detector underground (~ 100 m) and by 
only utilizing upward going slepton decays. The background from upward going neutrino 
interactions, giving rise to upward going muons, is eliminated by considering only upward 
going decay products originating from inside a fiducial volume defined to be within the 
outer RPC and TGC layers. 

Our analysis shows that definite evidence of the existence of the decays of trapped 
sleptons would be available in the first year of running with a received luminosity of 10 
fb~^. A measurement of the lifetime to 4% is possible after 3 years at nominal luminosity 
and running conditions and a received luminosity of 120 fb~^. An observation of a trapped 
slepton signal in this manner would presumably justify the more extensive and expensive 
water traps and instrumented resevoirs discussed elsewhere. 
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2. Slepton Properties 



Current collider bounds from the null searches for long lived tracks at LEP-II |20l require 
the slepton mass (?ti^-) to be greater than 99 GeV. Cosmology also has something to say 
about the slepton mass. In the gravitino super WIMP dark matter scenario, assuming that 
superWIMPS provide most, if not all, of the cold non-baryonic dark matter with density 
in the range 0.094 < QoMh'^ < 0.124, the mass of the slepton is constrained to the range 
m,- ~ 700-1000 GeV. Sleptons this heavy will be difficult to study at the LHC. 

However, if gravitinos are produced during reheating, then for the reheating tempera- 
tures preferred for leptogenesis as discussed above ( ~ 10^ GeV), gravitinos could form all 
of the non-baryonic dark matter for slepton masses as low as 120 GeV Alternatively, 
if the neutralino (x) is slightly heavier than the slepton, then it will freeze out at the same 
time as the slepton and then decay to the slepton, thus combining its relic density with that 
of the slepton. This permits much lower slepton masses to produce the correct gravitino 
dark matter density. For example, in minimal supergravity (mSUGRA) with Aq = 0, tan/3 
= 10, and /i > 0, the required relic density may be achieved in the region tlspXlsp at 

When the gravitino is not the LSP, 
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this is in the excluded stau LSP region. 

Assuming the lepton mass is negligible, the width of the slepton decay I 



IG is: 
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(2.1) 



As the gravitino decays gravitationally, the decay is determined only by the lepton mass, 
the gravitino mass and the Planck mass - no SUSY parameters enter. 

In the scenario considered in this paper where the gravitino is the LSP, the lower 
bound on Mi 12 (300 GeV) is determined by the requirement of a stau NLSP. The number 
of trapped staus generated rapidly diminishes as increases and an upper bound on 

Ml/2 ( = 800 Gev/c2) can be defined to be where only a few staus may be trapped in 
ATLAS per year. 

Knowledge of the slepton range in matter is a key aspect of this analysis. Charged 
particles passing through matter lose energy by emitting radiation and by ionizing atoms. 
In this analysis, we assume the slepton loses energy in a similar way to a very heavy 
muon of the same mass, where ionization losses dominate over energy loss by radiation. In 
this case, the GEANT3 simulation package |18] was used to estimate the range. Another 
author |l^ directly estimated the energy loss of the slepton per gcm~^ using the Bethe 
Bloch formula. The results of the GEANT simulation are in excellent agreement with the 
direct calculation for the range of a stau. 



3. The Model of the Detector and its Surroundings 

The ATLAS detector and experiment caverns, and the ATLAS rock overburden and un- 
derburden are all modelled in GEANT 3.21 [18|. Our detector model, shown in Figure ^, is 
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based closely on the ATLAS detector |23], as described in the ATLAS Technical Design Re- 
ports [^]. This model includes the beam pipe and shielding, the inner detector (the pixel 
detector, the semiconductor tracker (SCT) and the transition radiation tracker(TRT)), 
the calorimetry (the electromagnetic calorimeter (ECAL), the hadronic tile calorimeter 
(HCAL), the liquid argon endcap calorimeter (LAr) and the forward calorimeter (FCAL)), 
the magnet system (the barrel solenoid and the barrel and endcap toroids) and the muon 
system (the monitored drift tube chambers (MDT), the cathode strip chambers (CSC), the 
resistive plate chambers (RFC) and the thin gap chambers (TGC)). 

All volumes are filled with a material that is a smeared average of the materials found 
in the corresponding actual detector volumes, where these averages correspond to the 
radiation lengths of each of the ATLAS components. An ideal toroidal field is used to model 
the magnetic field inside the barrel and endcap toroid models, and the barrel solenoid has 
a uniform field of 2 T. 

The muon system, particularly the RPCs and TGCs, has the most direct impact on 
this study of trapped stau decay. The RPCs and TGCs are dedicated trigger layers and, in 
combination with the precision of the MDTs, can allow the muon system to stand on its own 
as a muon detector. The RPCs lie in the barrel section of the detector and are separated 
into three layers. The central layer is referred to as the pivot layer, where hits in this layer 
are projected outward to the inner and outer RPC layers in coincidence windows (cones 
about the direction of the incoming muon that are dependent on the muon's transverse ^ 
momentum) to search for matching hits. The outer layer is used only to search for high 
transverse momentum muons. The TGCs are split into five layers in each of the ATLAS 
endcap sections. All but one of these layers is in the form of a doublet, containing two 



^The transverse direction is orthogonal to the direction of the beam. 
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Figure 1: The range in lead of a stau with mass 219 GeV c ^ as a function of energy, estimated 
using a GEANT3 simulation of the stau range. 
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carbon dioxide-n-pentane gas volumes that are each layered between graphite cathodes 
and separated by paper honeycomb. The other TGC layer is a triplet, consisting of three 
gas volumes. The outermost doublets are the endcap pivot layers, where hits there are 
projected inward to the other TGC layers and matching hits are searched for in the same 
fashion as the RPCs. 



4. The Monte Carlo Simulation 

An R-parity conserving mSUGRA SUSY breaking model is assumed, with the common 
scalar mass nio and the common soft breaking parameter Aq taken to be GeVc, the ratio 
of the Higgs vacuum expectation values tan/3 is taken to be 10, the sign of the higgsino mass 
term /x taken to be greater than 0, and the common gaugino mass varied between 

300 GeV and 800 GeV in increments of 100 GeV. This one-dimensional parameter set lies 
in a region of the mSUGRA parameter space that requires the gravitino to be the LSP, 
otherwise it would be an excluded stau LSP region [19|. The lower bound on M1/2 



chosen to have a stau NLSP, and the upper bound is chosen so that the SUSY cross-section 
is high enough to have enough staus created each year to make this analysis useful. Table 
|l] lists the SUSY cross-section that corresponds to each value Mi/2- 



Ml/2 (GeV) 


300 


400 


500 


600 


700 


800 


crsuSY (pb) 


20.2 


4.87 


1.51 


0.542 


0.219 


0.096 



Table 1: The SUSY cross-sections that correspond to the values of AI1/2 considered. 




Figure 2: GEANT model of an LHC detector based closely on ATLAS. 
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To study trapped stau decay, the proton-proton collisions of the LHC are simulated 
using HERWIG 6.507 |5| interfaced with ISAJET 7.71 |6| in order to obtain the SUSY 



mass spectrum. We consider sparton, gaugino and/or slepton production generated using 
IPROC=3000 [p7| where the three following processes are implemented: 



• 2-parton — )• 2-sparton processes. All QCD sparton, i.e. squark and gluino, pair 
production processes are implemented; 

• 2-parton — )• 2-gaugino or gaugino+sparton processes. All gaugino, i.e. chargino 
and neutralino, pair production processes and gaugino-sparton associated production 
processes are implemented; and, 

• 2-parton — )• 2-slepton processes. All Drell-Yan slepton production processes are im- 
plemented. 

The stau lifetimes that we consider are 7, 30 and 90 days. Each stau is given a random 
creation time throughout a year during those periods when the LHC beam is on, and a 
random decay time according to its lifetime spectrum. 

After being assigned a lifetime, creation time and decay time, the staus are passed 
through the GEANT ATLAS model. GEANT tracks the staus until the stau momentum 
falls below 0.001 GeV/c^, at which point it is assumed the stau has become trapped. It 
is also assumed that all staus will stop before decaying ^. Also, because the stau is heavy 
and has an electromagnetic but not a colour charge, GEANT treats it like an extremely 
massive muon as it travels through matter, losing most of its energy through ionization 
effects. 

Once the stau has stopped in ATLAS, it is decayed by GEANT, and the resulting 
muons and pions down the decay chain are tracked as they pass through ATLAS. At that 
point, a model trigger decision is applied that mimics the decision made by the cosmic ray 
muon trigger. Based on |^8[, it is assumed that the trigger in the barrel region has an 83% 
efficiency, and the trigger in the endcap region has a 95% efficiency. 



5. Triggering on Trapped Slepton Decay 

Once the LHC attains its design luminosity of 10^^ cm~^s~^, ATLAS will become subject 
to more than 10® events per second. Because of their high cross-sections, QCD processes 
will dominate these events and any rare signals, like SUSY processes, will be swamped by 
this background. However, ATLAS utilizes a three-tier trigger system to extract these rare 
events and successively reduce the data rate to a final recordable rate of approximately 200 
Hz. 

The low- level trigger, level 1 (LVLl), is hardware-based and consists of custom elec- 
tronics operating at the LHC bunch crossing frequecy of 40 MHz with a latency time of 
just 2.5 fis. Up to 256 distinct trigger types that rely on information from the calorimeters 

^This is valid because the lifetimes considered here are much greater than the time required for the stau 
to travel even 1 km away from ATLAS through the rock overburden or underburden. 
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and muon detectors may be programmed into the LVLl trigger menus. The level 2 (LVL2) 
trigger and event filter (EF) make up the software-based high-level trigger that runs on 
dedicated processor farms. 

The ATLAS cosmic ray trigger consists of a set of dedicated LVL2 cosmic ray muon 
trigger algorithms that rely heavily on the LVLl muon trigger architecture. In the barrel 
region, the LVLl cosmic ray muon trigger |29] requires coincidences in both the rj and cj) 



directions of the inner two RFC layers. If the RFC pivot layer registers a hit, then the 
trigger searches for a coincidence in the inner RFC layer across all sectors of the inner 
RFC layer that are electrically connected. This loosens the pointing requirement that is 
inherently built into most LVL2 trigger algorithms, and corresponds approximately to the 
muon pointing to within 2 m of the interaction point. The timing of the cosmic rays is 
determined using the MDT chambers. 

The physical cross-section of the TGC endcaps is very small with respect to the cosmic 
ray muons that arrive from above essentially parallel to the TGC surface. Consequently, 
to boost the number of triggers that can occur, dummy hits are introduced into the large 
wheels to force a trigger to fire when only a single layer of the TGC triplet is actually hit. 
Since the layers of the TGC triplet are so closely spaced, any particle originating from the 
interaction point will likely register a hit in all three layers, with a track that may be traced 
back to the interaction point and that has a matching track in the inner detector. Cosmic 
ray muons will generally not point to the interaction point, so they are easily distinguished 
as cosmic rays. 

According to ATLAS, then, the signal of the trapped staus will look very similar to 
that of cosmic ray muons. The stau decay products will not generally originate from the 
interaction point, nor will they fit nicely into the bunch crossing periods of the LHC. As 
such, the ATLAS cosmic ray muon trigger is a natural choice to search for the signal of 
the trapped stau decay. Because the cosmic ray trigger can perform precision timing of 
the muons, it will also work for decay products originating both from within and outside 
of ATLAS for staus trapped either within ATLAS or in the rock surrounding the main 
ATLAS experiment cavern. 

6. Backgrounds 

In order to eliminate the overwhelming background produced by the proton-proton col- 
lisions, as well as the background from beam-halo and beam-gas rates [^^, we propose 
to search for the signal of trapped stau decay only when the LHC beam is off. The 
only remaining background sources to the trapped stau signal arise from cosmic rays and 
upward-going neutrino-induced muons. Neutrino-induced upward-going staus may also be 



a possiblity |31[| . 



ATLAS is bombarded by cosmic ray muons at a rate of kHz, which would likely com- 
pletely shadow the trapped stau signal. This large flux can be immediately dealt with 
by discarding all downward-going tracks, which is well within the capacity of the ATLAS 
trigger. While this would cut the number of signal events roughly in half, the signal to 
background ratio will be significantly enhanced by only considering upward-going events. 
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Ml/2 (GeV) 


7 days 




30 days 


90 days 




oUU 


713.6 




14.0 


709.5 




34.8 


693.8 




zl.o 


400 


673.1 




72.4 


705.1 




23.0 


718.9 




20.3 


500 


722.8 




29.0 


706.3 




27.1 


702.2 




15.2 


600 


713.8 




37.1 


702.5 




27.8 


685.5 




19.9 


700 


703.5 




32.0 


708.8 




29.2 


716.4 




8.2 


800 


706.3 




33.7 


702.6 




20.3 


693.0 




13.5 



Table 2: Acceptance (xlO^) of the ATLAS-like detector, defined as the total number of upward- 
going stau decay products observed by the model trigger relative to the total number of Monte-Carlo 
staus produced, from staus that are trapped within the fiducial volume. 

The background from upward-going muons created in interactions in the rock beneath 
ATLAS can also be ehminated. The "ATLAS-like" cosmic ray trigger considered here 
requires that the particle must pass through the inner two RPC layers and within 2m of 
the interaction point (in the barrel region), or it must pass through the innermost TGC 
layer (in the endcap). In order to rule out muons entering ATLAS from below, we restrict 
the trapping region for the staus comprising our signal to be inside the outer RPC layer 
(within a radius of 1027.5 cm) and inside the inner TGC layer (inside a z half-length of 
1291 cm). We also require the particles to be upward going. Looking only at particles 
trapped within the above-defined fiducial volume, we can be sure that we are not seeing 
anything but signal. 

7. Slepton Lifetime Determination 

The ratio of the number of stau decays seen by the model trigger to the total number of 
staus created gives the geometric detector acceptance, a. Table |2| gives a for upward-going 
staus trapped within the fiducial volume as a function of both M1/2 ^-nd lifetime, where a 
total of 10^ staus were produced at each point and a is an average over six different trials. 

The number of candidate stau decays that ATLAS will be sensitive to during a period 
when the beam is off will be proportional to the geometric acceptance a multiplied by 
the stau production rate. Thus, the number of candidate stau decays ATLAS will detect 
during a shutdown period as a function of elapsed time well be: 

Nhits = aNlje-'/\ (7.1) 

where r is the stau lifetime, t is the amount of time since the beam shut off and j is the 
number of trapped staus at beam shut off (see the Appendix for the complete definition). 
To increase the size of the data set, a sum over multiple years can be taken, leading to 

A^h^i. = a |^f;iV/^.e-*/^j , (7.2) 

where n is the number of years considered. 
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This proportionality may be expressed as 



AT = Ae-^^ (7.3) 



n 



A = aY,Kj (7.4) 

i=i 

B = \. (7.5) 

As such, it is possible to model-independently determine the lifetime of the stau using 
equation |7.3| to fit the distribution of the total number of candidate stau decay products 
seen as a function of elapsed time during a beam shut off period. The SUSY cross-section 



is embedded in the value of A (see equation |7.4|) , and may thus be determined, although 



the depedence of A on a makes this model-dependent. 

Two scenarios are considered in this analysis. The first assumes the LHC has run for 
only one year at low luminosity (an integrated luminosity of 10 fb~^). The second scenario 
assumes that three years worth of data have been collected, with two years of LHC running 
at low luminosity (as above) and one year of LHC running at high luminosity (an integrated 
luminosity of 100 fb~^). An analysis completed after only one year of LHC running could be 
used to demonstrate whether a decay signal may be seen. Combining three years worth of 
data would allow for a more accurate determination of the lifetime of the decaying particle, 
while still within a reasonable amount of time to see results. 



The LHC operating schedule is taken from [35[ |36|, which results in 200 days of beam 
running, segmented by 3 day beam shut off periods and a large 144 day beam shut down 
period over the winter. We use this large shut down period^, where the assumed lifetimes 
of 7, 30 and 90 days are shorter than the total amount of time the beam is off. 

To determine the statistical spread in these values, this fit was performed on 400 
distinct data sets for each scenario, using lifetimes of 7, 30 and 90 days, over the entire 
range of Mx/2- We only attempted this fit when the expected number of events was greater 
than ten with the bin size set so that there were at least 4 events in the first bin. Figure ^ 
shows, as an example, a distribution of the number of stau decays versus the elapsed time 
during the shutdown period, summing over 3 years of data (2 years at low luminosity and 
1 year at high luminosity). The solid curve depicts the fit to the data that gives the stau 
lifetime as 29.9 it 1.5 days (model independent) and the SUSY cross-section as 18.9 ± 1.3 
pb. Figure Q shows the distribution of lifetimes obtained for 400 experiments, under the 
second scenario, for a true stau lifetime of 30 days at M1/2 = 300 GeV/c^. In Figure |5|, 
we see the corresponding distribution of number of events for the 400 experiments, and in 
Figure ^ the distribution of the SUSY cross-section values. 

Tables ^, ^ and ^ give the results for the both scenarios where the slepton lifetimes 
are 7 days, 30 days and 90 days, respectively. One can see for all scenarios that, as the 
value of Mii2 increases, the SUSY cross-section decreases and, consequently, so does the 
accuracy of the stau lifetime and SUSY cross-section determination using trapped staus. 



^If the lifetime of the stau is not long compared to the 3 day shut off periods, then those periods could 
be used to perform this analysis, as well. This possibility is not considered here. 
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Figure 3: Sum over 3 years of data (2 years 
low luminosity and 1 year high luminosity) 
of the number of stau decay products seen 
by the model trigger as a function of elapsed 
time during the shutdown period. The fit of 
the distribution results in a lifetime of 29.9 
± 1.5 days and a SUSY cross-section of 18.9 
± 1.3 pb. 



Figure 4: The lifetime distribution for 400 
experiments under the second scenario where 
the true stau lifetime is 30 days and the value 
of Ml/2 is 300 GeV/c^. 




400 420 44C 450 4S0 500 520 540 

Njnber of detected staus 



Figure 5: The distribution of the number of 
events for 400 experiments under the second 
scenario where the true stau lifetime is 30 
days and the value of M1/2 is 300 GeV/c^. 
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Figure 6: The SUSY cross-section distribu- 
tion for 400 experiments under the second 
scenario where the true stau lifetime is 30 
days and the value of M1/2 is 300 GeV/c^. 



For r = 30 days, near the cosmologically preferred value of the stau lifetime (37 days), 
our results show (see Table § ) that under Scenario 1 for Mi 12 ~ 300-400 GeV and SUSY 
cross-section (T5(75y=20.2-4.87 pb the stau lifetime and asusY determined agree with the 
actual values within a fractional (statistical) error on the measured value of ~ ib(25%-60%). 
Under this scenario we obtain ~ 40-10 signal events after 1 year running at low luminosity. 
For scenario 2, for ~ 300-600 GeV with (T5[/5y=20.2-0.54 pb we can obtain a lifetime 
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and cross-section estimate that agrees with the set values within a fractional statistical 
uncertainty on the measured value of of ~ =b(5%-60%) and ~ =b(7%-70%), respectively, 
with ~ 467-13 signal events detected. 

Under Scenario 1, for a relatively short stau lifetimes of 7 days, we only have sensitivity 
for Ml/2 ~ 300 GeV with cross-section crsc/sy =20. 2 pb. In this case, the stau lifetime 
and SUSY cross-section determined agrees with the actual value to within the fractional 
(statistical) uncertainty of ~ ±50% and ~ ±80%, respectively, on the measured value, with 
~ 10 signal events detected. We can do better under scenario 2 where, for ~ 300-500 
GeV and asuSY = 20.2-1.51 pb, we can determine the stau lifetime and asusY to within 
the fractional (statistical) uncertainty on the measured value of ~ ±(12%-65%). In this 
case the number of detected events varies between ~ 120 and ~ 10. 

The longest stau lifetime considered is 90 days. For Scenario 1 we have sensitivity for 
Ml/2 ~ 300-400 GeV/c^ and (^SUSY = 20.2-4.87 pb. Under this scenario, we can determine 
the stau lifetime (Tsusy to within the fractional (statistical) uncertainty on the measured 
value of ~ ±(23%-30%) and ~ ±(13%-24%), respectively, with the number of recorded 
signal events varying between ~ 83 and ~ 21. For scenario 2, we can measure the stau 
lifetime over the range M1/2 ~ 300-700 GeV/c^ and asuSY = 20.2-0.22 pb. In this case, we 
can pin down the stau liftime and (Tsusy to within the fractional (statistical) uncertainty 
on the measured value of ~ ±(7%-40%) and ~ ±(4%-48%), respectively, with the number 
of detected events varying between ~ 990 and ~ 12. 

8. Conclusion 

The sensitivity of the ATLAS-like detector to the decay of long-lived staus that become 
trapped in the detector was studied using a simplified model of the ATLAS detector. To 
completely eliminate the background from cosmic rays and upward going neutrino-induced 
muons, only upward-going decay products originating from inside a fiducial volume defined 
to be within the outer RFC and TGC layers were studied. In addition, decays were only 
recorded during beam-off periods to further reduce backgrounds from actual collisions. 

By fitting the number of candidate stau decay products detected by a model cosmic 
ray muon trigger during periods when the LHC beam is off as a function of the elapsed 
time during the beam shut off period, the model-independent stau lifetime and the model- 
dependent SUSY cross-section can be determined. 

The fitting procedure returns the stau lifetime and SUSY cross-section with up to 
approximately 5% statistical accuracy on the measured value. For almost all situations 
considered that allow a fit - those with ~ 10 detected events or more - the measured values 
of the slepton liftime and SUSY cross-section lie within ±lcr of the actual value, where a 
is the statistical error on the measured value. We found that for the two running scenarios 
considered, our analysis procedure has a good sensitivity to slepton lifetimes of ~ 30 days 
and a SUSY cross-section greater than approximately 0.50 pb. It is interesting to note that 
stau lifetimes (~ 30 days are cosmologically favoured by models propounding a gravitino 
LSP scenario. 
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Table 3: Fitted and expected values of stau lifetime, SUSY cross-section and number of events 
for a stau lifetime of 7 days at both luminosity scenarios. The unit of lifetime is days and that of 
cross-section is picobarns. 



Appendices 



If we set the luminosity for year j of running at the LHC to be Lj, then the stau production 
rate at ATLAS, assuming R-parity conservation, will be 

R, = 2^2^, (1) 

J- on 

where Ton is the total time over a year that the beam is on. The schedule of beam operation 
assumed is given in |35]|36|, which gives a total of 200 days of LHC beam operation in one 



year. The integrated luminosity of the LHC is assumed to be 10 fb^^ and 100 fb^^ for low 
and high luminosity years, respectively. 

Now, letting the number of undecayed staus remaining in or around ATLAS at beam 
turn-on be N° and the number remaining at beam turn-off be A^-^, the number of staus 
remaining as a function of time t becomes 

Noff = A^-^e"*/^ 

Non = N"e-'/^ + RT{l-e-"^), (2) 

where r is the stau lifetime. 
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Table 4: Fitted and expected values of stau lifetime, SUSY cross-section and number of events 
for a stau lifetime of 30 days at both luminosity scenarios. The unit of lifetime is days and that of 
cross-section is picobarns 

Using equation |2| and labelling the eight LHC operating periods with the index i gives 
Nil = 

N?j = Nl^_^e~'''^l\ j > 2 

= iV°j-e~*°"/^ + RjT{l - 6^*°"/^), (3) 

where tsd, toff and ton are the amounts of time during a shutdown, offline and operational 
period of the LHC ^. This set of equations can be used to determine the number of un- 
decayed staus remaining at the beginning and end of each operational period, resulting 
in 

JSf^. = iV°^.e~*(*°//+*°")/'^e*o///^ -1- 

i-l 

RjT{l - e-*°"/^) e-"(*°//+*°")/^ 

n=0 

i~l 7 

m=l ra=0 

*The shutdown period is a long 144 day period when the beam is shut off and the offline periods are 
short three day down times between beam operational periods. 
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Table 5: Fitted and expc;ctc(l values of stau lifetime, SUSY cross-section and number of events 
for a stau lifetime of 90 days at both luminosity scenarios. The unit of lifetime is days and that of 
cross-section is picobarns. 

Meanwhile, the number of staus that decay during shutdown period i of year j is 

^decays ^ ^/^j _ g-W,/r)^ ■ ^ 7] 

^decays ^ ^/^^j _ ^-t^a/ry (5) 

Thus, the number of decays ATLAS will sec will be the the total number of decays that 
occur during a period when the beam is off multiplied by the geometric detector acceptance 
a. 
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